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Abstract. I review recent calculations of the suppression of bottomonium states in heavy 
ion collisions. A non-relativistic potential is used which is complex valued. This allows one to 
extract the binding energies and decay widths of the ground and excited states of bottomonium 
as a function of the typical plasma particle momentum and momentum-space anisotropy. The 
decay widths determined arc used as input and integrated over space-time taking into account 
the dynamical evolution of the typical particle momentum and momentum-space anisotropy. 
The suppression of T(ls), T(2s), T(3s), Xbi, and X62 is obtained as a function of centrality, 
rapidity, and transverse momentum. The obtained results are compared with data from the 
STAR and CMS collaborations. 



1. Introduction 

In November of 2010 the Large Hadron Collider (LHC) at the European Organization for Nuclear 
Research (CERN) achieved its first relativistic heavy ion collisions of lead nuclei (Pb-Pb). The 
goal of such experiments is to recreate conditions which only existed in the early universe before 
the formation of hadrons. The LHC is searching for evidence of the creation of a primordial 
state of matter called the quark gluon plasma (QGP) which comprised the entire universe until 
approximately 10~^ seconds after the big bang when the temperature of the universe was on 
the order of 10^^ Kelvin. The LHC experiment is attempting to recreate such temperatures 
by colliding lead nuclei at i/s/v7v = 2.76 TeV/nucleon whereas the lower energy Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Laboratory has been colliding gold nuclei 
at -^/sjviv = 200 GcV/nucleon. In addition to simply crossing the threshold for the creation of 
a quark gluon plasma, the experiments are also making detailed measurements of the collision 
products in order to determine fundamental properties of the QGP. The resulting data are being 
compared to theoretical calculations based on quantum chromodynamics (QCD) which describes 
the interactions of quarks and gluons. 

At high temperatures one expects the emergence of Debye screening of the interaction between 
quarks and gluons. This leads to the dissolution of hadronic bound states [1]. A particularly 
interesting subset of hadronic states consists of those comprised of heavy quarks, since the 
spectrum of low lying states can be found using potential-based non-relativistic treatments. 
Based on such potential models there were early predictions [2,3] that J/i/j production would 
be suppressed in heavy ion collisions relative to the corresponding production in proton-proton 
collisions scaled by the number of nucleons participating in the collision. In recent years there 
have been important theoretical advances in the understanding of heavy quark states at finite 
temperature using analytic techniques [4-12] and lattice QCD [13-21]. Most recently, interest 
has shifted to bound states of bottom and anti-bottom quarks (bottomonium) for the following 
reasons 



(i) Bottom quarks (m& 2± 4.2 GeV) are more massive than charm quarks {rric ~ 1.3 GeV) and 
as a result the heavy quark effective theories underpinning phenomenological apphcations 
are on much surer footing. 

(ii) Due to their higher mass, the effects of initial state nuclear suppression are expected to be 
smaller than for the charmonium states [22] . 

(iii) The masses of bottomonium states {mr ~ 10 GeV) are much higher than the temperatures 
(T ^ 1 GcV) generated in relativistic heavy ion collisions. As a result, bottomonium 
production will be dominated by initial hard scatterings. 

(iv) Since bottom quarks and anti-quarks are relatively rare within the plasma, the probability 
for regeneration of bottomonium states through recombination is much smaller than for 
charm quarks. 

As a result one expects the bottomonium system to be a cleaner probe of the quark gluon 
plasma than the charmonium system for which the modeling has necessarily become quite 
involved. For this reason we will focus on the bottomonium states and only consider the thermal 
suppression of these states, ignoring initial state effects and any possible thermal generation or 
recombination. In this conference proceedings we will review recent theoretical calculations of 
bottomonium suppression at energies probed in relativistic heavy ion collisions at RHIC and 
LHC. We will present an overview of the important aspects of the calculation and refer the 
reader to Refs. [23, 24] for details. We will compare our predictions for inclusive T(ls) and 
T(2s) suppression with recent data to from the CMS and STAR collaborations. 

2. Theoretical methods 

In the last few years there have been important theoretical developments in the theoretical 
treatment of heavy quarkonium in the quark gluon plasma. First among these are the first- 
principles calculations of imaginary-valued contributions to the heavy quark potential. The 
first calculation of the leading-order perturbative imaginary part of the potential due to gluonic 
Landau damping was performed by Lainc ct al. [5]. Subsequently, an additional imaginary- 
valued contribution to the potential coming from singlet to octet transitions has also been 
computed using the effective field theory approach [9] . These imaginary- valued contributions to 
the potential are related to quarkonium decay processes in the plasma. The consequences of such 
imaginary parts on heavy quarkonium spectral functions [25,26], perturbative thermal widths 
[5,27], quarkonia at finite velocity [28], in a T- matrix approach [4,8,29-31], and in stochastic 
real-time dynamics [32] have recently been studied. 

Additionally, there have been significant advances in the dynamical models used to simulate 
plasma evolution. In particular, there has been a concerted effort to understand the effects 
of plasma momentum-space anisotropics generated by the rapid longitudinal expansion of 
the matter along the beamline direction. The resulting dynamical models are now able to 
describe the anisotropic hydrodynamical evolution using full (3+l)-dimensional simulations [33- 
38]. This is important because momentum-space anisotropics can have a significant impact 
on quarkonium suppression since in regions of high momentum-space anisotropy one expects 
reduced quarkonium binding [7,10-12,39]. In Refs. [23,24] the dynamical evolution of the 
anisotropic plasma was combined with the real and imaginary parts of the binding energy 
obtained using modern complex-valued potentials. 

2.1. Anisotropic potential model and binding energies 

Early on it was shown that ideal relativistic hydrodynamics is able to reproduce the soft collective 
flow of the matter and single particle spectra produced at RHIC [40-43]. Based on this, there was 
a concerted effort to develop a systematic framework for describing the soft collective motion. 



This effort resulted in a number of works dedicated to the development and application of 
relativistic viscous hydrodynamics to relativistic heavy ion collisions [44-62]. 

One of the weakness of the traditional viscous hydrodynamics approach is that it relies on 
an implicit assumption that the system is close to thermal equilibrium which implies that the 
system is also very close to being isotropic in momentum space. However, one finds during 
the application of these methods that this assumption breaks down at the earliest times after 
the initial impact of the two nuclei due to large momentum-space anisotropics in the px-PL 
plane which can persist for many fm/c [63]. In addition, one finds that near the transverse 
and longitudinal edges of the system these momentum-space anisotropics arc large at all 
times [35,36,63]. Similar conclusions have been obtained in the context of strongly coupled 
systems where it has been shown using the conjectured AdS/CFT correspondence one achieves 
viscous hydrodynamical behavior at times when the system still possesses large momentum- 
space anisotropies and that these anisotropics remain large throughout the evolution [64-70]. 
Based on these results one is motivated to apply a dynamical framework that can accommodate 
potentially large momentum-space anisotropies. 

In order to take into account plasma momentum-space anisotropy, the phase-space 
distribution of gluons in the local rest frame is assumed to be given by [7, 71-74] 

fit, X, p) = /i,o (v'p'+e(p-n)VPhard) , (1) 

where f[so is an isotropic distribution which in thermal equilibrium is given by a Bosc-Einstcin 
distribution, ^ is the momentum-space anisotropy parameter, and Phard is a momentum scale 
which specifies the typical momentum of the particles in the plasma and can be identified with 
the temperature in the limit of thermal isotropic (C = 0) equilibrium. The two parameters phard 
and ^ can, in general, depend on proper time and position; however, we do not indicate this 
explicitly for compactness of notation. 

In general one finds that the heavy quark potential has real and imaginary parts, V = 
?li[V]+i'^[V] . One can determine the real part of the heavy-quark potential in the non-relativistic 
limit from the Fourier transform of the 00-component of the static gluon propagator. In an 
anisotropic plasma with a distribution function given by Eq. (1) one finds [10], at leading order 
in the strong coupling constant, 

V{r,0 = -^e'^-^A^%u = 0,p,O, (2) 



^ ^ I (27r)3 (p2 + ml + m2)(p2 + m}) - mj ' ^ ' 

where g is the strong coupling constant and Cp = {N^ — \)/{2Nc) is the quadratic Casimir of 
the fundamental representation of SU{Nc). The mass scales ma-, m^, m^, and mg are listed in 
Ref. [10]. One can factorize the denominator of (3) by introducing 

2m% = M^± - 4(m^(m2 + m^) - mj) , (4) 

with = -I- -|- [71]. This allows us to write 

^ '^^ ^ y (27r)3 (p2 + m^)(p2 + m?.) 



In general one must evaluate (5) numerically. The integration can be reduced to a two- 
dimensional integral over a polar angle, 9, and the length of the three-momentum, p. However, 



there can be poles in the integration domain due to the fact that m?_ can be negative for 
certain polar angles and momenta [71]. These poles are first order and can be dealt with 
using a principle-part prescription. After evaluating the integral numerically, one finds that 
the resulting potential is very well described by a Debye-screened Coulomb potential with an 
anisotropic Debye mass /x 

nnr)] = -'-^e--^\ (6) 

where 

\mn) ' + H (3 + 0'' )[ + ) ' 

The coefficients a-c are fixed by (a) requiring that the small and large anisotropy limits of the 
analytic potential are reproduced and (b) fitting to the numerical results obtained by direct 
integration of Eq. (5) [24]. 

Here we will focus on a model in which the real part of the potential is obtained from internal 
energy of the system since models based on the free energy seem to be incapable of reproducing 
either the LHC or RHIC i?AA[T"]- The real part of the potential is given by [24] 

mV] = -"(l + //r)e-'^^ + — [l-e-'^n - f^^e-'^'' - ^> (8) 
r jJL f^Qf 

where a = 0.385 and a = 0.223 GeV^ [75] and the last term is a temperature- and spin- 
independent finite quark mass correction taken from Ref. [76]. In this expression n = Q{^, 9)mD 
is an anisotropic Debye mass where ^ is a rather complicated function (see Eq. (7) above) 
which depends on the degree of plasma momentum-space anisotropy, ^, and the angle of the 
line connecting the quark-antiquark pair with respect to the beamline direction, 9, and is 
the isotropic Dcbyc mass [24]. In the limit ^ — > one has = 1. In the figures in the results 
section, results obtained with Eq. (8) are often labeled as "Potential Model B" . 

The imaginary part of the potential S[V^] is obtained from a leading order perturbative 
calculation which was performed in the small anisotropy limit 

'^[V] = -asCFT{<l){r/mD) - C [Mr/nin, 9) + ^2{r/mDM} , (9) 

where 0, ^i, and %Ij2 can be expressed in terms of hypergeomctric functions [39]. After combining 
the real and imaginary parts of the potential, the 3d Schrodinger equation is solved numerically 
to obtain the real and imaginary parts of the binding energy as a function of ^ and Phard [12, 77]. 
The imaginary part of the binding energy is related to the width of the state 

^ r\-i 2SJ[£^bind(r,Xi,0] 5?[-Ebind(T,Xx,0] > 0, ..^x 

[r, x^, g - I 5R[£;bind(T, x^, <j)] < , ^'^^ 

where <; = arctanh(z/t) is the spatial rapidity. The value of 10 GeV in the second case is chosen 
to be large in order to quickly suppress states which are fully unbound (which is the case when 
the real part of the binding energy is negative). 



2.2. Dynamical Model 

The dynamical model used gives the spatio-temporal evolution of the typical transverse 
momentum of the plasma partons, Phard(7"7x), and the plasma momentum-space anisotropy, 
4(r, x), both of which are specified in the local rest frame of the plasma. The widths obtained 
from solution of the 3d Schrodinger equation are then integrated and exponentiated to compute 
the relative number of states remaining at a given proper time. This quantity is then averaged 
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Figure 1. Real and imaginary parts of the T(ls) (left) and T(2s) (right) binding energies as a 
function of the hard momentum scale, phard for different levels of momentum-space anisotropy 

over the transverse plane taking into account the local conditions in the plasma and weighting by 
the spatial probability distribution for bottomonium production which is given by the number 
of binary collisions computed in the Glauber model with a Woods-Saxon distribution for each 
nucleus. For the temporal integration the initial time is set by the formation time of the state 
in question. The resulting Raa is a function of the transverse momentum, pT, the rapidity 
and the nuclear impact parameter b. To compare to experimental results transverse momentum 
cuts are applied assuming a n{pT) = noE^'^ spectrum. In addition, any cuts on the rapidity due 
to detector acceptance and centrality are applied as needed. For details of the dynamical model 
and Raa computation we refer the reader to Ref. [24]. 

2.3. Initial conditions 

For the initial conditions we use a Woods-Saxon distribution for each nucleus and determine 
the transverse dependence of the initial temperature via the third root of the number of 
participants (wounded nucleons). In the spatial rapidity direction we have investigated two 
possible temperature profiles: (a) a broad plateau containing a boost-invariant central region 
with Gaussian limited-fragmentation at large rapidity [78] 

n(?) = no exp(-(|<?| - aat/2)V2fT,') ©(kl " ?fiat/2) , (11) 

where <jflat = 10 is the width of the central rapidity plateau, cr^ = 0.5 is the width of the limited 
fragmentation tails, and no is the number density at central rapidity [57]; and (b) a Gaussian 
motivated by low-energy fits to pion spectra 

n(?) = noexp(-?V2cr^) with cr^ = 0.64 • 8 In (^/i^/2mp) /3(1 - ) , (12) 

where Cg = l/\/3 is the sound velocity, nip = 0.938 GeV is the proton mass, and ^/sNN is 
the nucleon-nucleon center-of-mass energy [24,79]. The temperature distribution is given by 
T ~ n^/^. We note that Eq. (11) has the advantage that it has been tuned to successfully 
describe the rapidity dependence of the elliptic flow in LHC heavy ion collisions. 

3. Results and Conclusions 

In Fig. 1 we show the real and imaginary parts of the T(ls) and T(2s) binding energies as a 
function of the hard momentum scale, Phard; for C ^ {0, 1, 20}. Defining the disassociation scale 




Figure 2. RHIC and LHC suppression factor Raa for the T(ls), T(2s), T(3s), xwi and Xb2 
states as a function of the number of participants (left) and rapidity (right). In the top plots 
we used y^sjvTV = 200 GeV and in the bottom plots we used ^/~SNN = 2.76 TeV. In both cases, 
we assumed a shear viscosity to entropy density ratio of Airrj/S = 1, and implemented cuts of 
< < 20 GeV and (left) rapidity \y\ < 0.5 for RHIC and \y\ < 2.4 for LHC (right) centrahty 
0-100%. 



as the value of Phard at which the real and imaginary parts of the binding energy become equal, 
one finds isotropic dissociation scales of 593 MeV for the T(ls) and 228 MeV for the T(2s). In 
addition, as one can see from this figure, as the level of momentum space anisotropy increases, 
one finds that the disassociation scale increases. This means that in regions of the plasma where 
there is a high degree of momentum-space anisotropy, the states will be less suppressed. 

In Fig. 2 we show the predicted suppression factor Raa for the T(ls), T(2s), T(3s), Xbi, 
and Xb2 states as a function of the number of participants (left) and rapidity (right). In this 
figure we see clear signs of sequential suppression, with the higher excited states having stronger 
suppression than the ground state. However, we note that even for states that are "melted" at 
relatively low central temperatures, we still obtain a non- vanishing suppression factor for these 
states. This is due to the fact that near the edges, where the temperature is lower, one does 
not see suppression of the states. Upon performing the geometrical average, we see that a large 
fraction of the states produced can survive even when the central temperature of the plasma is 
above their naive dissociation temperature. 

In Fig. 3 (left) we plot i?yi^[T(ls -|- 2s -|- 3s)] and compare with experimental data from the 




Figure 3. RHIC T(ls + 2s + 3s) suppression factor (left) compared witli experimental data 
from the STAR Collaboration [80] . The three different lines correspond to different assumptions 
for the shear viscosity to entropy ratio Airrj/S G {1,2,3}. In both plots we used ^/sN^ = 200 
GeV and implemented cuts of < < 20 GeV and |y| < 0.5. On the right, the solid black 
line is the result obtained assuming wounded nucleon initial conditions and the dashed red line 
is the result obtained used a two component model with a = 0.145. 



STAR Collaboration [80]. As can be seen from this figure, the model does a reasonably good 
job of reproducing the existing STAR data for Raa[^{^s + 2s + 3s)]. From the left panel we can 
obtain an estimate of rj/S: 0.08 < rj/S < 0.24. In Fig. 3 (right) we show the results obtained 
for i?^^[T(ls + 2s + 3s)] at RHIC energies for two different types of initial conditions (Glauber 
and mixed Glauber plus binary collision scaling). For both lines shown in Fig. 3 (right) we have 
assumed Attij/S = 2. Because changing the initial condition type affects particle multiplicities 
we have adjusted the initial temperature at RHIC energies from 433 MeV to 461 MeV in order 
to keep the charged particle multiplicity fixed at dN^h/dy = 620. As can be seen from Fig. 3 
(right), for peripheral collisions there is a larger dependence on the choice of initial condition 
type, while for central collisions the result obtained is not much affected by the choice of initial 
condition. This is to be contrasted with the dependence on the assumed value of r]/S which 
affects the suppression at all centralities. 

When considering the suppression of the T(ls) and T(2s) states it is important to include 
the effect of feed-down from higher excited states. In pp collisions only approximately 51% 
of T(ls) states come from direct production and similarly for the T(2s). One can compute 
the inclusive suppression of a state using i?^]{[T(ns)] = Z^i e states /« -^i.Ayi where fi are the 
feed-down fractions and Ri^AA is the direct suppression of each state which decays into the 
T(ras) state being considered. Here we will use fi = {0.510,0.107,0.008,0.27,0.105} for the 
T(ls), T(2s), T(3s), Xbii and Xb2 feed-down to T(ls), respectively [82]. For the inclusive T(2s) 
production we use fi = {0.500,0.500} for the T(2s) and T(3s) states, respectively. For details 
of the computation of the direct Raa for each state see Ref. [24] . 

In Fig. 4 we compare to recent data on the inclusive T(ls) and T(2s) suppression available 
from the CMS collaboration [83]. For this figure we used a broad rapidity plateau as the 
initial density profile as specified in Eq. (11). The central temperatures were taken to be 
To = {520,504,494} MeV at tq = 0.3 fm/c for Airr^/S = {1,2,3}, respectively, in order to 
fix the final charged multiplicity to dNch/dy = 1400 in each case. As can been seen from this 
figure, the predictions agree reasonably well with the available data. The data seem to prefer the 




Figure 4. Predictions for the central rapidity inclusive (a) T(ls) and (b) T(2s) suppression 
including feed-down as a function of A'part along with recent data from the CMS collaboration. 
The three different lines correspond to different assumptions for the shear viscosity to entropy 
ratio Atttj/S £ {1,2,3}. 




Figure 5. LHC inclusive suppression factor Raa for the T(ls) including feed down effects as 
a function of rapidity and transverse momentum compared to experimental data are from the 
CMS Collaboration [81]. 



largest value of r]/S shown; however, there is a ±14%, ±21% Is, 2s global uncertainty reported 
by CMS, making it hard to draw firm conclusions. 

In Fig. 5 we show the inclusive suppression factor i?J^]|[T(ls)]. Comparing to the available 
CMS data [81] we can obtain an estimate for r]/S at LHC energies: 0.08 < ij/S < 0.24 which 
is the same range obtained from the STAR data obtained with gold-gold collisions at lower 
energies. As before, more precisely determining rj/S will require more data from the LHC which 
should be forthcoming in the near future. 

4. Conclusions 

In this proceedings contribution we have reviewed recent calculations of bottomonium 
suppression in RHIC and LHC heavy ion collisions. We presented comparisons of prior 



predictions with data from the CMS Collaboration for inclusive T(ls) and T(2s) suppression and 
data from the STAR Collaboration for inclusive T(ls + 2s + 3s) suppression. The underlying 
calculations employed a complex-valued potential which incorporates both screening and in- 
medium dissociation of the states under consideration. We solved for the resulting real and 
imaginary parts of the binding energy of the T(ls), T(2s), T(3s), Xbi^ ^'Hd Xb2 states. We then 
folded this information together with the real-time evolution of both the typical momentum 
of the plasma particles (phard) and their momentum-space anisotropy (^). We demonstrated 
that the resulting inclusive suppression of T states is in good agreement with available data. 
We note in closing that there are now very interesting developments concerning quarkonium in 
an anisotropic medium obtained employing the conjectured AdS/CFT correspondence [84-86]. 
These methods offer some hope to determine the temperature and momentum-space anisotropy 
dependence of the long range part of the potential from first principles, albeit in a theory that 
is not QCD. 
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